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1.  INTRODUCTION. 


Consider  the  problem  of  comparing  v  treatments  called  test  treatments  with  a 
standard  treatment  called  the  control  in  b  blocks  of  size  k  each.  The  test  treat¬ 
ments  will  be  labelled  l,2,...v  and  the  control  0.  The  model  is  additive  and 
homoscedastic,  that  is,  if  Yij|  is  the  observation  on  treatment  i  (0  <  i  <  v)  in 
block  j(l  <  j  <  b)  and  plot  1  (1  <  1  <  k),  then 

Yijl  =  V  +  ri  +  +  eijl> 

where  the  e^  are  assumed  to  be  uncorrelated  random  variables  with  mean  0  and 
a  common  variance  a2.  The  unknown  constants  ft,  rt  and  (3}  represent  the  general 
mean,  the  effect  of  treatment  i  and  the  effect  of  block  j  respectively.  Let 
D(v,b,k)  be  the  set  of  all  possible  experimental  designs.  The  primary  purpose  of 
the  experiment  is  to  draw  inferences  on  the  contrasts  (r0  —  Tj),  i  =  l,2,...,v.  If 
(rdo  -  i"di)  *s  the  best  linear  unbiased  estimator  of  (r0-r;)  under  a  design 

d  e  D(v,b,k),  then  we  want  to  choose  a  design  such  that  the  variances  of 

(rdo  -  rdi)  are  smallest  in  some  sense.  Formally,  we  want  to  choose  an  experimen¬ 

tal  design  from  D(v,b,k)  which  minimizes 

V 

E  var(rdo  "  *di) 
i=l 

as  d  varies  over  all  of  D(v,b,k).  A  design  which  attains  the  minimum  is  called  an 
A-optimal  design.  Throughout  this  paper,  we  will  assume  that  v  and  k  satisfy 

k  >  3  (1.1) 

v  >  k.  (1.2) 

Majumdar  and  Notz  (1983)  gave  a  method  for  finding  A-optimal  designs. 
Their  optimal  designs  can  basically  be  of  two  types.  Using  the  terminology  of 
Hedayat  and  Majumdar  (1984),  they  are:  rectangular  (or  R-)type,  in  which  every 
block  has  the  same  number  of  replications  of  the  control,  and  step  (or  S-)type,  in 
which  some  blocks  contain  the  control  t  times  and  the  others  t  +  1  times. 
Optimal  R-type  designs  were  studied  by  Hedayat  and  Majumdar  (1985).  Families 
of  such  designs,  particularly  when  each  block  has  one  replication  of  the  control, 
were  given  in  that  paper.  In  this  article,  we  intend  to  study  optimal  S-type 
designs.  S-type  designs  are  more  complicated  than  R-type  designs;  the  latter 
being  a  balanced  incomplete  block  (BIB)  design  in  the  test  treatments  augmented 
by  an  equal  number  of  controls  in  each  block,  but  the  former  does  not  have  such 
a  simple  characterization.  Consequently,  both  the  optimality  and  the  construc¬ 
tion  of  such  designs  are  more  involved. 


In  section  2  we  give  a  characterization  of  classes  D(v,b.k)  for  which  the 
Majumdar  and  Notz  (1983)  method  gives  an  optimal  S-type  design.  Hedayat  and 
Majumdar  (1985)  contains  a  similar  result  for  R-type  designs.  The  results  of  sec¬ 
tion  2  are  applied  in  section  3  to  obtain  some  families  of  optimal  S-type  designs. 
We  also  suggest  ways  of  obtaining  efficient  designs  in  classes  D(v,b,k)  where 
optimal  designs  are  not  readily  available. 


2.  THE  NATURE  OF  OPTIMAL  STEP  TYPE  DESIGNS. 

We  start  this  section  by  giving  some  definitions  and  known  results.  For 

d  e  D(v.b,k),  ndii  will  denote  the  number  of  times  treatment  i  (0  <  i  <  v)  occurs 
J  b 

in  block  j  (1  <  j  <  b),  and  Xdip  =  ^n^n^  (0  <  i,p  <  v).  The  following 
definition  is  due  to  Bechhofer  and  Tamhane  (1981). 

Definition  2.1.  d  is  a  Balanced  Treatment  Incomplete  Block  (BTIB)  design  if 
^doi  =  '  '  ~  ^dOv 

and 

^d!2  =  ^dl3  =  ■  ‘  '  —  ^d,v-I,v- 

A  special  type  of  BTIB  designs  will  be  of  particular  interest  to  us. 

Definition  2.2.  For  integers  t  €  {0,1,  •  •  •  ,  k  -  1},  s  6  {0,1,  •  •  ,  b  -  1},  d  is  a 

BTIB(v,b,k  ;  t,s)  if  it  is  a  BTIB  design  with  the  additional  property  that 

ndij  e  {0,1},  i  =  1,2,  •  •  •  ,  v;  j  =  1,2,  ■  •  ,  b, 

nd01  ~  =  Dd0s  =  t  +  1, 


ndo,s+ 1  —  ■  •  •  —  ndob  —  t. 

The  layout  of  a  BTIB(v,b,k  ;  t,s)  design  can  be  pictured  as  follows,  with  columns 
of  the  array  denoting  the  blocks: 


Here  di  and  d2  are  components  of  the  design  which  involve  the  test  treatments 
only.  If  s  =  0  then  d!  is  empty  and  d2  is  necessarily  a  BIB  design.  In  this  case 
the  BTIB(v,b,k  ,  t,0)  is  called  rectangular  (R-)type.  Otherwise,  when  s  >  0,  the 
BTIB(v,b,k  ;  t.s)  is  called  step  (S-)type. 


a  =  (v  -  l)2,  c  =  bvk(k  -  1),  p  =  v(k  -  1)  +  k, 

A  =  { ( x , z ) :  x  =  0,  •  ■  ,  [k/2)  -  1;  z  =  0,1,  ■  •  ,b.  with  z  >  0  when  x 


Here  [  -  )  is  the  largest  integer  function.  Let, 
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t 

I 

l 

g(x,z)  =  a,  {c  -  p(bx  +  z)  -t-  (bx2  -t-  2xz  +  z)} 

j  4-  1  / { k ( bx  -+-  z)  -  (bx2  +  2xz  +  z)}. 

Observe  that  g(x,0)  =  g(x  -  l,b),  for  all  x  0.  The  following  theorem  is  due  to 
Majumdar  and  Notz  (1983). 

Theorem  2.3,  A  BTIB(vrb,k  ;  t,s)  is  A-optimal  in  D(v,b,k)  if 
g(t,s)  =  min  g(x,z). 

(x,i)eA 

oome  examples  of  optimal  R-  and  S-  type  designs  are  available  in  Majumdar 
and  Notz  (1983)  and  in  the  catalogs  of  optimal  designs  given  in  Hedayat  and 
Majumdar  (1984).  It  may  be  noted  that  the  designs  given  in  Theorem  2.3  are 
also  Maximum  Variance  (MV-)  optimal,  that  is  they  minimize 

max  var(fd0  -  rdi) 

1  <  i  <  v 

among  all  deD(v,b,k).  More  general  results  on  MV-optimal  designs  are  now 
available  in  Jacroux  (1984).  Hedayat  and  Majumdar  (1985)  investigated  the 
nature  of  g(x,z)  and  characterized  classes  D(v,b,k)  for  which  the  g(x,z)  is  minim¬ 
ized  at  (t,0)  for  some  t. 

Let  L  be  the  set  of  integers  r  with  0  <  r  <  |k/2]b.  Then  the  function  g  in 
Theorem  2.3  can  also  be  expressed  as  a  function  G  on  L  as  follows: 

G(r0)  =  g(x,z),  where  r0  =  bx  +  z  with  x  =  |r0/b)  and 

z  =  r0-  b ( r o ,  b|. 

Thus  a  BTIB(v,b,k  ;  t,s)  is  A-optimal  if  G(r0)  is  minimized  at  r0*  =  bt  +  s.  Here 
r0  is  the  number  of  observations  allocated  to  the  control.  For  specific  values  of 
v.b  and  k,  the  minimization  of  G  over  L  can  be  performed  on  a  computer 
without  difficulty.  However,  if  one  wants  to  derive  some  general  results  such  as 
the  optimality  of  a  family  of  designs,  (for  instance,  Theorem  3.1  in  the  present 
paper)  then  it  would  be  necessary  to  study  some  theoretical  properties  of  the 
function  G  (or  g).  The  following  result,  first  conjectured  in  Ture  (1982),  provides 
an  important  technical  tool  for  investigating  the  A-optimality  of  BTIB  designs: 

Theorem  2.4.  If  G(r0*)  =  minG(r0),  then  G  is  decreasing  on  (r0  £  L:  r0  <  r0*} 

r^L 

and  increasing  on  {r0  €  L:  r0  >  r0*}. 

A  Proof  of  Theorem  2.4  can  be  found  at  the  end  of  this  section. 

Let  Aq  =  {(0,z) :  z  =  1,  •  •  •  ,  b}  and  Aj  =  {(i,z) :  z  =  0,1,  •  •  •  ,  b}  for 
i  =  1,2,  •  •  •  ,  (k/2|  -  1.  Then,  by  Theorem  2.4,  one  of  the  following  two  situa¬ 
tions  must  prevail,  for  some  t  €  {0,1,  •  •  •  ,  (k/2)}: 

a)  g(x,z)  is  a  decreasing  function  of  z  on  each  of  A^A,,  •  •  ■  ,  At_1  and,  if 
t  <  |k/2|,  g(x,z)  is  increasing  on  each  of  At,  •  •  ■  ,  A|k/2j_|.  In  this  case,  the 


minimum  of  g(x,z)  over  A  is  achieved  at  x  =  t  and  z  =  0,  and  an  R-type  design, 
BTIB(v,b.k  ;  t,0)  is  optimal.  This  was  the  case  studied  by  Hedayat  and  Nlajum- 
dar  (1985). 

b)  g(x,z)  has  a  local  minimum  at  (t,s)  on  At  for  some  s  €  {1.  ■  •  ■  .  b  -  1 },  g(x,z)  is 
decreasing  on  each  A;,  i  <  t,  if  any,  and  g(x,z)  is  increasing  on  each  A;,  i  >  t.  if 
any.  In  this  case,  an  S-type  design,  BTIB(v,b,k  ;  t,s)  is  optimal. 

An  immediate  consequence  of  Theorem  2.4  is  that  a  local  minimum  of  G  is 

also  a  global  minimum.  Thus  G(r0*)  =  minG(r0)  if  and  only  if 

r<jeL 

G(r0*)  <  min{G(r0*  -  1),  G(r0*  +  1)}.  In  other  words,  we  have 

Corollary  2.5.  For  some  t  €  {0,1,  •  •  •  ,  (k/2|  -  1}  and  s  €  { 1,  ■  ••  ,  b  -  1 }, 
g(t,s)  =  min  g(x,z)  if  and  only  if 

(x,j)€A 

g(t,s)  <  min{g(t,s  -  1),  g(t,s  +  1)}.  (2.1) 

Consequently,  a  BTIB(v,b,k  ;  t,s)  is  A-optimal  in  D(v,b,k)  whenever  (2.1)  holds. 

Condition  (2.1)  substantially  reduces  the  computation  one  has  to  do  in  verify¬ 
ing  the  optimality  of  a  BTIB  design.  In  section  3,  it  will  be  used  to  show  the 
optimality  of  a  family  of  S-type  BTIB  designs. 

Theorem  2.4  is  also  useful  for  simplifying  the  search  for  the  optimal 
r0*  =  bt  -+•  s,  which,  as  will  be  demonstrated  in  the  next  section,  is  useful  for  con¬ 
structing  efficient  designs  when  optimal  designs  cannot  be  found  among  BTIB 
designs.  One  simple  method  for  finding  the  optimal  (t,s)  could  be  the  sequence: 

(i)  Determine  x0  =  min{x:  g(x,b  -  1)  <  g(x,b),  x  =  0,1,  •  •  •  ,  |k/'2|  -  1}. 

(ii)  Determine  z0  such  that  g(x0,z0)  =  min  g(x0,z).  A  method  of  determining 

i=0,l,...,b 

z0  could  be  to  consider  g(x0,z)  as  a  function  of  a  real-valued  z  over  the  interval 
[0,b|,  and  using  techniques  of  Calculus  to  find  its  minimum.  Suppose  the 
minimum  is  at  a  point  p  €  |0,b).  If  p  is  an  integer,  then  z0  =  p.  If  p  is  not  an 
integer,  then  z0  =  Jp]  if  g(x0,[p))  <  g(x0,|p)  +  1),  z0  =  [p]  +  1  if 
g(x0,|p|)  >  g(x0,|p)  +  1),  and  any  one  of  [p]  and  |p)  +  1  can  be  chosen  as  z0  if 
g(xo>|p])  =  g(x0,(p!  +  1).  In  the  last  situation  a  choice  between  [p|  and  [p|  +  1 
should  be  made  with  an  eye  to  the  existence  of  an  optimal  BTIB  design. 

(iii)  If  z0  =  0  then  t  =  x0,  s  =  0.  If  z0  €  {1,  •  •  ■  ,  b  -  1}  then  t  =  x0,  s  =  z0.  If 
z0  =  b  then  t  =  x0  +  1,  s  =  0. 

Remark.  It  is  easily  seen  that 
g(x,z)  =  b'g(x,z), 

where  z  =  b_,z  and  g  does  not  depend  on  b.  Thus,  treating  z  as  a  real  variable 
on  (0,1),  one  can  now  minimize  g(x,z)  over  x  e  {0,1,  ■  •  •  ,  |k/2|  -  1}  and  z  e  [0,1) 
to  obtain  optimal  x0,  z0.  The  solutions  xn  and  i#  dfi  depend  op  1?,  If  z0b  is 
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an  integer  then  the  optimal  values  are  t  =  xQ  and  s  =  z0b.  Otherwise,  the 
minimum  of  g(x,z)  over  A  is  attained  at  x  =  x0  and  z  =  |z0b]  or  |z0bj  -  1.  This 
fact  makes  it  possible  to  tabulate  (x0,z0)  as  a  function  of  k  and  v.  One  such 
table  is  available  in  Ture  (1985). 

.As  an  example,  consider  v  =  6,  k  =  5  and  b  =  18.  Here  x0  =  1  and  the 
optimal  value  of  z  is  z0  =  .321,  but  z0b  =  5.778  is  not  an  integer.  Since 
g|l,6)  <  g(  1 .5),  the  optimal  s  =  6.  A  BTIB(6,18,5  ;  1,6)  is  given  in  Hedayat  and 
Majumdar  (1984)  (design  S4  in  its  appendix). 

We  finish  this  section  by  proving  Theorem  2.4.  For  this  purpose  we  need  two 
lemmas.  Lemma  2.6  was  first  proved  in  Ture  (1982);  also  see  Ture  (1985)  and 
Hedayat  and  Majumdar  (1985). 

Lemma  2.6.  The  function 
z  —  g(x,z) 

is,  for  a  fixed  x  6  {0,1,  •  •  •  ,(k  -  2 ) / 2 } ,  a  decreasing  function  on  |0,b|,  or  an 
increasing  function,  or  there  is  a  z0  €  (0,b)  such  that  the  function  is  decreasing  on 
(0,z0|  and  increasing  on  |z0,b). 

Lemma  2.7.  (i)  If  x  €  {0,1,  •  •  •  ,  (k/2)  -  2}  then  g(x,b)  >  g(x,b  - 1)  implies 

g(x  +  1,1)  >  g(x  +  1,0). 

(ii)  If  x  €  {1,2,...,  fk/2|  -  1}  then  g(x,I)  <  g(x,0)  implies 
g(x  -  1 , b )  <  g(x  -  l,b  -  1). 

Proof.  (i)  It  follows  from  p.762  of  Hedayat  and  Majumdar  (1985)  that 
g( x , b )  >  g(x,b  -  1)  if  and  only  if 

b{k  -  (x  +  1 ) } [ a{ x  +  l)2(p  -  2x  -  1)  —  (k  -  2x  —  l){v(k  -  1)  -  (x  +  l)}2) 

>  v(k  -  2x  -  l)(p  -  2x  -  l){k  -  l  +  (v  -  2)(x  +  1)}. 

In  order  that  the  above  inequality  holds,  the  left  hand  side  should  be  positive 
since  the  right  hand  side  is.  Hence  if  g(x,b)  >  g(x,b  -  1),  then 

a(x  +  l)2(p  -  2x  -  1)  >  (k  -  2x  -  l){v(k  -  1)  -  (x  +  l)}2.  (2.2) 

Again,  from  p.761  of  Hedayat  and  Majumdar  (1985),  g(x  +  1,1)  >  g(x  +  1,0)  if 
and  only  if 

b{k  -  (x  +  l)}[(k  -  2x  -  3){v(k  -  1)  -  (x  +  l)}2  -  a(x  +  l)2(p  -  2x  -  3)j 

<  v(k  -  2x  -  3)(p  -  2x  -  3){k  -  1  +  (v  -  2)(x  +  1 )}. 

This  inequality  is  certainly  valid  if  the  left  hand  side  is  negative,  since  the  right 
hand  side  is  positive;  that  is,  if 

a(x  +  l)2(p  -  2x  -  3)  >  (k  -  2x  -  3){v(k  -  1)  -  (x  +  1  )}2. 

It  suffices  to  show  that  (2.2)  implies  (2.3).  Using  (2.2)  we  see 


(2.3) 


a(x  -  l)2(p  -  2x  -  3)  =  a(x  +  1  )2(p  -  2x  -  1)  -  2a(x  -t-  l)2 

>  (k  -  2x  -  1 ){ v( k  -  1)  -  (x  +  l)}2  -  2a(x  ■+■  l)2 

=  (k  -  2x  -  3){v{k  -  1)  -  (x  +  l)}2  +  2{v(k  -  1)  -  (x  l)}2  -  2a(x  -  l)2. 
So,  it  suffices  to  show  that 

{v(k  -  1)  -  (x  +  l)}2  >  a(x  +  l)2, 

that  is  {v(k  -  1)  +  (v  -  2 )( x  +  l)}{v(k  -  1)  -  v(x  +  l)}  >  0,  which  is  true  since 
x  <  (k  2)  -  2. 

The  proof  of  (ii)  is  completely  similar  and  involves  the  following  steps: 
g( x ,  1 )  <  g( x ,0 )  is  equivalent  to 

b(k  -  x)((k  -  2x  -  1 ) { v(k  -  1)  -  x}2  -  ax2(p  -  2x  -  1)] 

>  v(k  -  2x  -  l)(p  -  2x  -  l){k  -  1  +  (v  -  2)x}. 

Hence,  (k  -  2x  -  1  ){v(k  -  1)  -  x}2  >  ax2(p  -  2x  -  1),  so  that 
(k  -  2x  -h  1 ) { v ( k  -  1)  -  x}2  >  ax2(p  -  2x  +  1).  Hence, 

b(k  -  x)(ax2(p  -  2x  +  1)  -  (k  -  2x  +  1 ){ v(k  -  1)  -  x}2j 

<  v(k  -  2x  +  l)(p  -  2x  +  l){k  -  1  +  (v  -  2)x}, 

that  is  g(x  -  l,b)  <  g(x  -  l,b  -  1). 

Now  we  are  ready  to  prove  Theorem  2.4. 

Proof  of  Theorem  2.4: 

If  the  minimum  of  g(x,z)  over  A  is  attained  at  s(  ie  (t,z)  with  z  =  0,  then  the 
result  follows  from  Lemma  2.4  of  Hedayat  and  Majumdar  (1985).  On  the  other 
hand, suppose 

g(t,s)  =  min  g(x,z), 

(x,j)€A 

where  s  {0,b } .  Then  since 

g(t,s)  <  g(t,s  +  1)  <  *  <  g(t,b  -  1)  <  g(t,b), 

we  obtain  from  Lemma  2.7  that 

g(t  +  1,0)  <  g(t  +  1,1). 

Now,  using  lemma  2.6, 

g(t  +  1,0)  <  g(t  +  1,1)  <  •  •  •  <  g(t  +  l,b). 

Continuing  this  chain  of  arguments,  we  see  that 


g(t,s)  =  min  g(x,z) 

*=t . (k/2|-l, 

a=0,l,...,b 

Similar  arguments  starting  from  g(t,s)  <  g(t,s  -  1)  show  that 


-  9  - 


3.  APPLICATIONS. 

An  important  application  of  Corollary  2.5  is  finding  families  of  optimal  step 
type  designs.  We  give  some  families  in  Theorem  3.1.  We  shall  also  discuss  situa¬ 
tions  where  Corollary  2.5  cannot  be  applied,  and  suggest  ways  of  finding  efficient 
designs. 

Theorem  3.1.  Let  o  >  3  be  a  prime  or  a  power  of  a  prime.  Then,  for  any  posi¬ 
tive  integer  7,  there  exists  a  BTIB(a2  -  1,  -t-  2)(a2  -  1),  a;  0,  7(0  +  l)(a2  -  1)) 

which  is  A-optimal  in  D(a2  -  1,  7(0  4-  2)(a2-l),a). 

Proof:  The  A-optimality  follows  from  Corollary  2.5  upon  verifying  condition 
(2.1). 

We  will  establish  the  existence  for  7  =  1;  for  larger  7  we  take  the  union  of  7 
copies  of  this  design.  We  denote  by  BIB(v,b,r,k,X)  a  BIB  design  with  the  parame¬ 
ters  v,b,r,k  and  X;  the  symbols  enjoying  their  standard  interpretation. 

Since  a  is  a  prime  or  a  prime  power,  there  exists  a 
BIB(a2,  a(a  +  l),o  +  l,o,l),  that  is  an  Euclidean  plane.  There  is  also  a 

BIB(o2  -  1,  (a  +  1 )( or2  -  2),  a2  -  2,  a  -  1,  a  -  2).  To  see  this,  start  with  an 

Euclidean  plane  based  on  a2  treatments  and  delete  one  treatment  from  it.  The 
remaining  design  has  a  +  1  blocks  of  size  a  -  1  and  a2  -  1  blocks  of  size  a.  Take 
o-2  copies  of  each  block  of  size  a  -  1,  and  replace  each  block  of  size  o  by  o 
blocks  of  size  a  -  1,  the  a  blocks  being  the  a  subsets  of  size  a  -  1  of  the  original 
block  of  size  a.  These  (a  -  2)(a  +  1)  +  o(a2-  1)  =  (a  +  l)(a2-  2)  blocks  give 
the  desired  BIB  design. 

To  get  a  BTIB(a2  -  1,  (a  +  2)(a2  -  1),  a;  0,  (a  4  1 )( o2  -  1 )),  augment  one 
replication  of  the  control  to  each  block  of  the  constructed 
BIB(o2  -  1,  (a  +  l)(a2  -  2),  a2  -  2,  a  -  1,  a  -  2)  and  add  to  this  the  blocks  of  a 

BIB(a2,  a(a  +  1),  a  +  1,  o,l),  in  which  the  control  appears  as  one  of  the  treat¬ 

ments.  Hence  the  theorem. 

Stufken  (1986)  has  more  families  of  optimal  step  type  designs.  There  again, 
condition  (2.1)  is  crucial  in  proving  the  optimality. 

The  problem  of  the  existence  of  a  BTIB(v,b,k;  t,s)  imposes  restrictions  on  the 
scope  of  Corollary  2.5.  As  in  BIB  designs,  the  parameters  should  satisfy  some 
necessary  conditions  (Hedayat  and  Majumdar  1984,  p.365).  Even  when  these  are 
satisfied,  a  BTIB(v,b,k;  t,s)  may  not  exist.  In  case  the  parameters  v,b,k  cannot 
accommodate  Corollary  2.5,  one  may  try  to  find  a  design  in  D(v,b,k)  which  is 
highly  efficient,  even  if  it  is  not  known  to  be  A-optimal.  It  is  not  difficult  to  see 
from  Theorem  2.3  that  for  any  d  e  D(v,b,k), 

V 

<r~2  Evar(rd0-  rdi)  >  vkg(t,s) 


(3.1) 


The  left  hand  side  of  this  expression  is  the  value  of  the  A-criterion  for  a  design  d. 
We  denote  this  by  Ad.  Clearly  the  ratio  vkg(t,s)/Ad  gives  a  measure  of  the 
efficiency  of  d. 

The  lower  bound  in  (3.1)  is  achieved  by  a  BTIB(v,b,k;  t,s).  Before  discussing 
specific  methods  of  identifying  efficient  designs,  let  us  determine  how  close  an 
arbitrary  design  in  D(v,b,k)  can  come  to  the  lower  bound  vkg(t,s).  For 
d  e  D(v,b,k),  let  us  defin? 

mdi  =  |{j  :  ndoj  =  l}j, 

that  is,  md|  is  the  number  of  blocks  in  which  the  control  appears  1  times, 
1  =  0,1,  •  •  •  ,  k  -  1.  Clearly 
k-l 

^lmdi  rdo, 

v  k-l 

-  Xdoi  =  rdo(k  -  1)  -  s  1(1  -  l)mdl, 

i=  l  1=1 

v  v 

and  if  ndii  €  {0,1}  for  all  i  =  1,2,  •  •  •  ,  v,  j  =  1,2,  •  •  •  ,  b,  then  ED  Adi  +  E  Adoi 

‘  i,p=  1  i— 

i<  p 

1  1  k-l  b 

=  — bk(k  -  1) - E  1(1  -  l)mdl  where  rdo  —  E  ndoj,  the  number  of  replications 

m  **  1 5=2 1  J  “  1 

of  the  control  in  d. 

Lemma  3.2.  Let  m^m,,  •  •  •  ,  mk_,  be  fixed  numbers  satisfying 
m0  +  mt  +  •  ■  •  +  mk_1  =  b.  Consider  the  subset  of  D(v,b,k)  in  which  md|  =  mlf 
1  =  1,2,  •  •  •  ,  k  -  1.  For  these  designs  denote  rdo  =  r0  and  Xdip  =  Aip, 
i,p  =  0,1,  ■  •  ■  ,  v.  Define  the  quantities 

X0  =  S\i/v,  X,  =  2E  E  Aip/v(v  -  1). 

i=i  i,p= i 


Then,  for  each  d  €  D(v,b,k), 

Ad  =  <r2  E  var(rdo  -  rdi)  >  vk(A0  +  Xi)/{X0(X0  +  vA!)}.  (3.2) 

i=  1 

Proof.  Follows  from  Lemma  2.2  of  Majumdar  and  Notz  (1983). 

The  condition  (3.2)  gives  a  sharper  lower  bound  to  Ad  than  (3.1).  It  can  be 
used  to  eliminate  configurations  (m^m],  •  *  •  ,  mk_i)  for  which  the  lower  bound  in 
(3.2)  is  not  very  close  to  vkg(t,s). 

To  find  an  efficient  design  when  Corollary  2.5  cannot  be  applied,  we  first 
determine  the  optimal  r0*  =  bt  +  s  by  the  method  described  in  section  2.  If 
there  did  exist  a  d  €  D(v,b,k)  which  was  a  BTIB(v,b,k;  t,s)  then  this  d  would 
have  been  A-optimal.  Since  there  is  no  such  design,  we  shall  look  at  values  r0 


close  to  r0*  and  various  configurations  of  (mQ.m,,  •  •  •  ,  mk_,).  For  each  of  these 
we  compute  lower  bounds  given  by  (3.2),  and  consider  those  which  are  not  much 
larger  than  vkg(t,s),  the  lower  bound  given  by  (3.1).  Next  we  consider  designs 
having  these  (m^m,,  •  •  •  ,  mk_,)  configurations  which  are  very  close'  to  BTLB 
designs.  By  this  we  mean  Xdoi's(l  <  i  <  v)  are  as  close  as  possible  and 
Xdip's(l  <  i,p  <  v,  i  7^  p)  are  as  close  as  possible.  We  choose  the  most  efficient 
design  among  all  these  designs.  It  is  expected  that  this  design  would  be  highly 
efficient,  in  general.  In  particular  we  expect  two  types  of  designs  to  perform  very 
well.  One  is  a  BTIB  design  with  rdo  close  to  r0*;  the  other  is  a  design  ‘closest’  to 
a  BTIB  design,  with  rdo  equal  to  r0*.  We  now  give  two  illustrative  examples. 

Example  3.3.  Let  v  =  5,  k  =  4  and  b  =  7.  Then  t  =  1  and  s  =  0,  r0*  =  7, 
vkg(  1,0)  =  2.04. 

Table  3.1 


(m0lm,,m2lm3) 

(X0,Xj) 

lower  bound  from  (3.2) 

(I, 6, 0,0) 

(3.6,  2.4) 

2.137 

(1,5, 1,0) 

(3.8,  2.2) 

2.134 

(0,7, 0,0) 

(4.2,  2.1) 

2.041 

(1,4, 2,0) 

(4,  2) 

2.143 

(0,6, 1,0) 

(4.4,  1.9) 

2.060 

(0,6, 0,1) 

(4.2,  1.8) 

2.165 

(0, 5,2,0) 

(4.6,  1.7) 

2.091 

Note  that  whenever  X0  and  X,  are  integers,  the  corresponding  design  achieving 
the  lower  bound  is  a  BTIB  design.  Consider  five  designs  d1,d2,d3,d4,d5,  with 

V 

columns  as  blocks,  and  the  corresponding  values  of  Ad  =  a  2  £  var(rdo  -  rdi): 

i=l 


0000002 
.  _  0111113  . 

ai  —  2223344  ’ 

3454555 


=  2.156 


0000000 

1111122 

2233433 

4545545 


-  2.058 


0000000 
0111222 
I  1  3  3  4  3  3  4 


=  2.067 


d4  = 


0000001 
0  0  1  1  2  2  2 
1  3  3  4  3  4  3 
2455554 


=  2.143  (BTIB  design) 


dc  = 


0000000 
0111123  » 

0222334  ’  d* 
>5  5  3  4  4  4  5 


=  2.207 


Clearly  d2  is  the  best  among  all  these  designs.  In  fact  Adj  is  even  smaller  than 
the  lower  bounds  given  in  table  3.1  for  all  other  configurations  of  (m0,m1,mo,m3). 
d2  is  not  a  BTIB  design  but  rdj0  =  r0*.  It  is  highly  efficient  since 

vkg(t,s)/Adj  =  .992. 

Note  that  the  BTIB  design  d4  has  a  high  efficiency  also,  since 
vkg(t,s)/Ad<  =  .952. 


L_3J 

Li  V  —  6,  k  =  5, 

b  =  7.  Here 

t  =  1  and  s  =  2, 

=  2.204. 

Table  3.2 

ro 

(m0,m,,m2,m3,m4) 

(X0lX,) 

lower  bound  from  (3.2) 

7 

(0,7, 0,0, 0,0) 

(4.67,  2.8) 

2.236 

8 

(0,6,1, 0,0,0) 

(5,  2.6) 

2.214 

9 

(0,5, 2, 0,0.0) 

(5.33,  2.4) 

2.204 

10 

(0,4, 3, 0,0,0) 

(5.67,  2.2.) 

2.207 

11 

(0,3, 4, 0,0,0) 

(6,  2) 

o  ooo 

12 

(0,2, 5, 0,0,0) 

(6.33,  18) 

2.249 

0000000 

0111122 

1223333 

2445445 

3566566 


=  2.243 


0000000 

0011112 

1222243 

3533454 

4656665 


2.227 
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(0  0  0  0  0  0  0 
0001112  . 

d3  =  1  1  2  2  3  4  3  ,  j\j  =  2.223 

2  3  3  4  5  5  4 
'5  4  6  6  6  6  5 


0000000 

0000113 

d4  =  1  1  2  2  2  2  4  ,  Aj  =  2.222  (BTIB  design) 

3434355 

5665466 

Clearly  the  BTIB  design  d4  is  the  best.  Note  that  rd<0  is  close,  though  not  equal 
to  r0*.  This  design  is  highly  efficient,  in  fact, 

vkg(t,s)/Ad<  =  .992. 

It  may  be  suggested  that  one  way  of  finding  efficient  designs  is  to  look  for  an  A- 
optimal  design  among  all  BTIB  designs  only.  Hedayat  and  Majumdar  (1984) 
shows  that  even  though  these  designs  perform  very  well,  in  general,  there  can  be 
instances  where  it  is  quite  poor.  This  is  because  the  stringent  combinatorial  con¬ 
ditions  of  BTIB  designs  may  enforce  a  very  inefficient  choice  of  r0.  For  instance, 
they  show  that  in  D(10,80,2)  it  is  possible  to  achieve  at  least  a  24%  improvement 
over  the  best  BTIB  design.  In  this  case,  the  optimal  (x0,z0)  =  (0,  0.49),  so 
bz0  =  80x0.49  =  39.2.  The  best  BTIB  design  has  r0  =  80,  which  is  too  far 
away  from  39.2.  We  recommend  to  use  the  right  r0  (or  approximately  so)  and 
then  construct  a  design  combinatorially  close  to  a  BTIB  design. 
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